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Abstract We discuss the implications of the recent X- 
ray and TeV 7-ray observations of the PSR B1259-63 
system (a young rotation powered pulsar orbiting a Be 
star) for the theoretical models of interaction of pulsar 
and stellar winds. We show that previously considered 
models have problems to account for the observed be- 
haviour of the system. We develop a model in which the 
broad band emission from the binary system is produced 
in result of collisions of GeV-TeV energy protons acceler- 
ated by the pulsar wind and interacting with the stellar 
disk. In this model the high energy 7-rays are produced 
in the decays of secondary neutral pions, while radio and 
X-ray emission are synchrotron and inverse Compton 
emission produced by low-energy (< 100 MeV) electrons 

1 from the decays of secondary charged 7r ± mesons. This 
model can explain not only the observed energy spectra, 

. but also the correlations between TeV, X-ray and radio 
emission components. 

Keywords pulsars : individual: PSR B1259-63 • 
X-rays: binaries • X-rays: individual: PSR B1259-63 
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1 Introduction 

PSR B1259-63 is a ^48 ms radio pulsar in a highly ec- 
centric (e~0.87), 3.4 year orbit with a Be star SS 2883 
P. The pulsar crosses the Be star disc twice per orbit, 
just prior to and just after periastron. Unpulsed radio, 
X-ray and 7-ray emission observed from the binary sys- 
tem are produced due to the collision of pulsar wind 
with the wind of Be star. Observations of the temporal 
and spectral evolution of the non-thermal emission from 
the system provide a unique opportunity to probe the 
physics of the pulsar winds (PW) which is, in spite of 

ISDC, Ch. d'Ecogia 16, 1290 Versoix, Switzerland 
E-mail: andrii. neronov@obs.unige.ch 



the wealth of observational phenomena, and a 40-year 
old observation history, remains a matter of debate. 

The interaction of the PW with the wind from the 
companion star, is responsible for the formation of a 
" compactified" pulsar wind nebula (PWN) with the the 
size about the binary separation distance (typically, on 
AU-scale). Compact size, large matter density and the 
presence of a strong source (companion star) which illu- 
minates the nebula make the physical properties of the 
compact PWN significantly different from the ones of 
their larger scale cousins. 

We present the results of the last observation cam- 
paign of the PSR B1259-63 system during the 2004 
pulsar periastron passage and their applications for the 
theoretical modelling of the source. We show that most 
of the observed properties of the system in radio-to- 
TeV band can be naturally explained within a model 
of proton-loaded pulsar wind. 



2 Multi-wavelength observations of the system 
during 2004 periastron passage. 

The upper panel of Fig. shows the X-ray lightcurve of 
the system [2] together with the TeV |3] and radio @] 
lightcurves. For comparison we show also the data from 
archival X-ray [SBS] and radio [7||5] observations. Rapid 
growth of the X-ray flux found in XMM-Newton obser- 
vations of 2004 is correlated with the rapid growth of the 
unpulsed radio emission from the system. The growth of 
radio and X-ray flux at these phases can be attributed 
to the pulsar entering the Be star disk. 

Unfortunately, TeV observations start somewhat later 
and it is not possible to see whether the TeV flux grows 
during the pre-periastron disk crossing. However, simple 
geometrical argument tells that the orbital phase 9 at 
which the pulsar should enter the disk for the second time 
should be shifted by 180° relative to the first entrance. 
From Fig.^one can infer that the first pre-periastron en- 
trance falls roughly between the phases 70° < 6 < 110°. 
Thus, the pulsar has to enter the disk again between the 
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Fig. 1 Comparison between the X-ray (top), TeV (middle) 
and radio (bottom) lightcurves. XMM-Newton observations 
are marked with triangles, BeppoSAX ones with circles, and 
ASCA ones with squares. Data for four different periastron 
passages, are shown with different colors: red (1994), green 
(1997), black (2000) and blue (2004). Bottom X axis shows 
the orbital phase, 9, top X axis shows days from periastron, 



phases 250° < 8 < 290°. Surprisingly, one can clearly see 
from the middle panel of Fig. ^ that the TeV flux grows 
in this phase interval. To test the conjecture that the TeV 
flux grows during the second entrance to the disk we su- 
perimpose the pre-periastron X-ray and TeV lightcurves 
over the post-periastron lightcurves by shifting the phase 
of the post-periastron data points by —180°. The result 
is shown in Fig. [21 One can see that in such representa- 
tion the rise and decrease of both X-ray and TeV flux 
from the system can be well fitted with a gaussian curve 
F(8) ~ exp (-(0 - 8 Q f/{2A8l)). We find that the best 
fit is achieved with the parameter choice 8o ~ 109.1°, 
Ado — 18.5° (the coordinate A8 along the X-axis of Fig. 
13 is, in fact A8 = 8 — 8q). The position of the Gaussian 
with the above parameters is shown schematically in Fig. 
n]by a shaded area. Denser and lighter shadings in Fig. 
2] correspond to the one and two widths of the Gaussian. 
Depending on the physical mechanisms of the X-ray and 
TeV emission, the inferred width of the Gaussian Ado 
gives either an estimate of the "thickness" of the Be star 
disk, or of the characteristic cooling time of the high- 
energy particles injected at the phase of the disk passage 
(see Section |5J). 

The graphical representation of the evolution of the 
X-ray photon index (the spectrum is well fit by an ab- 
sorbed powerlaw model) along the orbit is given in Figure 
|3J The most remarkable feature of the spectral evolution 
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Fig. 2 The X-ray (top) and TeV (bottom) flux as a function 
of the relative phase A9 = 9 — 9o (see text for the definition 
of 9o). The curves show a fit with a gaussian of the half- width 
A9 = 18.5°. 
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Fig. 3 Evolution of the X-ray photon index r p h over the 
orbital phase 9. Radio flux and spectral index evolution from 
the 1997 periastron passage |7j are shown in black. 



of the system is the hardening of the X-ray spectrum 
close to the moment when pulsar enters the Be star disk 
at the phase 8 ~ 8q — 2A8o — 70°. One can see that the 
decrease of the photon index r is simultaneous with the 
onset of the rapid growth of the X-ray flux. Similar hard- 
ening of the spectrum down to the photon index r p h ~ 1 
(or, equivalently, down to the spectral index a ~ 0) at 
the moment of disk entrance is observed in the radio data 
shown in Fig. in black 7 . To the best of our knowl- 
edge, neither the hardening of the X-ray spectrum, nor 
strong correlation between the radio and X-ray flux and 
spectral index variations was predicted in any of existing 
models of X-ray emission from the system. 
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3 Implications for theoretical models. 

Termination of the pulsar wind in the stellar disk leads 
to the acceleration of ultrarelativistic electrons and sub- 
sequent X-ray and gamma-ray emission. Within a sim- 
ple geometrical picture [Sj , the maximum of X-ray (syn- 
chrotron) emission is expected during the phases of the 
disk passage, while the maximum of TeV 7-ray(IC) emis- 
sion is expected at the moment of the periastron passage. 
However, surprisingly, the TeV 7-ray lightcurve (middle 
panel of Fig. ^) has a local minimum in the periastron. 
The observed X-ray - TeV and radio - X-ray correla- 
tions do not agree with the early theoretical predictions 
pillUj . To find the range of possible theoretical models 
which can explain the data it is useful first to make ba- 
sic qualitative estimates of different time scales present 
in the system. 



3.1 Characteristic cooling times. 

Electrons. One of the main differences between the syn- 
chrotron and IC mechanisms of X-ray emission is the 
difference in the cooling time scales. The cooling time of 
the TeV electrons which can produce synchrotron emis- 
sion at the energies 65 ~ 1 — 10 keV in the magnetic field 
B is 

f s ~6x 10 2 [S/O.lGf 372 [es/10 keVp 1/2 s. (1) 

The spectrum of optically thin synchrotron emission from 
the cooled electron population has the photon index r p h > 
1.5. Any hardening of the X-ray spectrum down to the 
values r p h < 1.5 (e.g. due to the increased injection of 
electrons at higher energies) would be "washed out" by 
the synchrotron cooling at the at the 10 2 — 10 3 s time 
scale. To the contrary, the typical IC cooling time in X- 
rays is 

tic (t) — 6 x 10 5 [i?/10 13 cm] 2 [e /c /10 keVp 1/2 s. (2) 

(we have assumed that the seed photons for the IC scat- 
tering come from the companion star of luminosity ~ 
10 38 erg/s and temperature T ~ 2 x 10 4 K; the subscript 
"(T)" indicates that the estimate applies for the Thomp- 
son regime). Estimating the size of emission region to be 
about the binary separation distance, R ~ 10 13 cm, one 
can find that electrons emitting IC radiation at 1 keV 
cool at the day time scales. Observation of the gradual 
evolution of the X-ray photon index down to r p h ~ 1.2 
and then back to r p h > 1.5 on the time scale of several 
days during the first entrance to the disk (see Fig. |3J) is 
consistent with the IC, rather than synchrotron model 
of X-ray emission. 

In principle, it is possible that X-ray and TeV 7-ray 
emission from the system are, respectively, low- and high- 
energy tails of the IC spectrum. Substituting naively the 
energy of TeV photons tic ~ 1 TeV into Eq. © one 
finds that the time scale of the spectral variability at TeV 



energies should be very short. However, at TeV energies 
the IC scattering proceeds in the Klein-Nishina regime 
and the cooling time in this regime grows with energy, 

tic (KN) =i 8.5 x 10 3 [tic /I TeV] ' 7 [i?/10 13 cm] 2 s (3) 

The minimum of the IC cooling time, ~ 10 3 s, is reached 
at roughly at the energy of transition between Thompson 
and Klein-Nishina regimes, E e ~ ejc ~ 10 — 100 GeV. 

Pulsar wind electrons are able to escape from the 
region of the dense photon background along the contact 
surface of pulsar and stellar wind. If the two winds do 
not mix, the pulsar wind flows along the contact surface 
with the speed wpw ~ 10 10 cm/s and escapes beyond the 
binary separation distance over the time scale 

t esc ~ R/vpw ~ 10 3 [i?/10 13 cm] s (4) 

This time scale is essentially shorter than the IC cool- 
ing time both in X-ray and in the TeV energy bands. 
Electrons escaping from the innermost region of pulsar 
wind/stellar wind interaction fill the larger extended re- 
gion (a " compactified" PWN of the size i?pwN of about 
several binary separation distances) and can loose their 
energy via IC emission at longer time scales in the less 
dense photon background produced by the Be star. The 
escape time from the compact PWN can be naively es- 
timated assuming diffusion in the weak PWN magnetic 
field. E.g. taking the diffusion coefficient D equal to the 
Bohm diffusion coefficient at E e ~ 1 TeV and depending 
on the energy as D ~ E~ a (a = 1 for the case of Bohm 
diffusion) one finds 

i PWN ~ 10 4 [5/0.1G] [E e /1 TeVp Q [i? PW N/T0 13 cm] 2 s(5) 

During the periods of the pulsar passage through the 
dense equatorial disk of Be star (typical density of the 
slow equatorial stellar wind at the location of the pulsar 
is ^disk ~ 10 10 — 10 11 cm -3 ), bremsstrahlung and ionisa- 
tion energy losses can compete with the IC loss. Indeed, 
the energy independent bremsstrahlung loss time, 

hrems = 10 4 [n/10 n cm" 3 ]^ s (6) 

is comparable to the IC loss time for the TeV electrons 
© and is shorter than the IC loss time for the X-ray 
emitting electrons J2J). Thus, during the short period of 
escape from the dense equatorial disk (the escape time 
is given by Eq. (J2J) as much as 

L~jL e ~ t esc /t brems ~ 10% [?idi S k/10 n cm -3 ] (7) 

of the power in relativistic electrons, L e , can be con- 
verted into the (bremsstrahlung) 7-ray luminosity, L 1 . 

Below the electron energy E e ~ 350 MeV the Coulomb 
energy loss, for which the cooling time scale is given by 

tcoui ^ 3 x 10 3 [n d isk/10 n cm -3 ]" 1 [£ e /100 MeV] s(8) 

dominates over the bremsstrahlung loss. During the pe- 
riods of the disk passage, essentially 100% of the power 
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output in electrons with energies below the "Coulomb 
break" 

Ecoui =s 30 [n disk /10 n cm- 3 ] MeV (9) 

(estimated from the condition that the Coulomb loss 
time is equal to the escape time, tcoui ~ t esc ) will be 
channeled into the heating of the disk, rather than on 
emission from the system. As a result, only electrons 
with energies above Ecoui can be injected into the com- 
pactified PWN. 

3.2 Protons. 

GeV-TeV energy protons can loose their energy only in 
interactions with the protons from the stellar wind. The 
enhancement of the pp interaction rate is expected dur- 
ing the pulsar passage through the dense equatorial disk 
of Be star. The pp interaction time 

tpp ~ 1.6 x 10 4 [n/10 11 cm" 3 ] _1 s (10) 

is comparable to the electron bremsstrahlung loss time 
(JBJ . Following the same way of reasoning as in the case of 
bremsstrahlung, one can find that as much as 10% of the 
power L p contained in the PW protons can be channeled 
in the secondary particles (7-rays, neutrinos, electrons, 
positrons) produced in pp interactions. The 7r° decay 7- 
rays carry away about 1/3 of the power output in pp in- 
teractions. Thus, the "7-ray efficiency" of pp interactions 
is somewhat lower than the efficiency of bremsstrahlung, 

L^/Lp ~ 0.3t esc /tpp ~ 3% [n disk /10 n cm- 3 ] (11) 

However, if the PW is proton-dominated, the luminosity 
of the 7-ray emission from pp interactions can exceed the 
bremsstrahlung luminosity. 



4 IC model of X-ray to TeV emission. 

Taking into account that the seed photons for the IC 
scattering have energies of about 10 eV (assuming the 
temperature of Be star T ~ 2 x 10 4 K), one can find 
that the IC emission from electrons of the energy E e 
peaks at 

ejc ^ 4 [E e /W MeV] 2 keV (12) 

The energy of the upscattered photons becomes approx- 
imately equal to the energy of electrons at 

e/c,(T^KN) ^ 30 GeV (13) 

(the transition to the Klein-Nishina regime) . If the spec- 
trum of electrons is a simple powerlaw with the spectral 
index p e (dN e /dE ~ E~ p "), the IC spectrum below and 
above the Thompson - Klein-Nishina break is, respec- 
tively, dN^/dE ~ S^+i)/ 2 and dN y /dE ~ E-P-+ 1 InE. 
The IC emission in the 10-100 GeV energy band is char- 
acterized by one more spectral feature. Namely, the IC 
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Fig. 4 IC model fit for the X-ray to TeV spectrum of the 
source. Radio emission is synchrotron from electrons which 
produce X-ray IC flux. 

cooling time of the 10-100 GeV electrons is comparable 
to the escape time from the compact region with a dense 
photon background. Estimating the energy of the cool- 
ing break in the IC emission spectrum from the condition 
t esc ~ tie, (T) one finds e/c, cool — 4 [i?/10 13 cm] GeV. 
Taking into account the coincidence of the cooling break 
energy with the energy of transition to the Klein-Nishina 
regime one can not expect to detect the conventional 
steepening of the IC spectrum above the cooling break 
because of the reduced efficiency of the IC scattering in 
the Klein-Nishina regime. One more complication of the 
detailed calculation of the IC emission spectrum in the 
GeV- TeV energy band is that in order to explain the 
observed behaviour of the TeV lightcurve during the pe- 
riastron passage within the IC model one has to assume 
that either additional non-radiative cooling mechanism 
dominates electron energy loss close to the periastron, or 
a cut-off in the electron spectrum at sub- TeV energies is 
present m] . The combined effect of the above mentioned 
difficulties makes the detailed predictions for the IC spec- 
trum in the GeV- TeV band quite uncertain and we do 
not attempt the detailed fit of the observed spectrum in 
this band. Instead we concentrate of the attempt to fit 
the general shape of the spectral energy distribution in 
the X-ray to TeV 7-rayband within the IC model. 

Fig.0]shows an example of the fit the the spectrum of 
PSR B1259-63 in IC model for X-ray to TeV emission. 
One can see that EGRET upper limit on the flux from 
the system requires the presence of a break in the IC 
spectrum at the energies E ~ 1 MeV. In the model fit 
shown in the Figure, the electron spectrum below the 
break at E e — 100 MeV has the spectral index p e = 2, 
while above the break the spectrum steepens to p e + 1 = 
3. It is clear that the overall shape of the IC spectrum in 
the keV to TeV energy band agrees well with the data. 

The energy of the break in the electron spectrum 
(~ 100 MeV) is close to the energy of the Coulomb break 
given by Eq. As it was discussed above, electrons 
with energies below Ecoui loose all their energy via the 
severe Coulomb loss before they are able to escape from 
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the dense equatorial disk of Be star to the less dense 
PWN. As a result, regardless of the initial injection spec- 
trum of electrons from the P W, the spectrum of electrons 
injected in the compact PWN has a low-energy cut-off 
at the energy ~ Ecoui ■ The IC cooling of electrons in the 
PWN leads to the formation of the characteristic power- 
law tail of electron distribution below Ecoui with p e = 2. 
The electron spectrum above the energy Ecoui (assumed 
to be a powerlaw with the spectral index p e = 3 in the 
model fit of Fig. 0J) is determined by the balance of ac- 
celeration and energy losses in the pulsar/stellar wind 
shock region. 

Electrons responsible for the X-ray IC emission pro- 
duce synchrotron radiation in radio band at the charac- 
teristic frequency 

e s ~ 1.5 [5/0.1 G] [E e /30 MeV] 2 GHz (14) 

The ratio of the synchrotron to IC luminosity is given 
by the ratio of the energy densities of the magnetic field 
and radiation, 

Ls/Lic = 2 x 1(T 4 [5/0.1 G] 2 [i? PWN /10 13 cm] 2 (15) 

The radio luminosity of the system is some 4 orders of 
magnitude lower than the X-ray luminosity. This im- 
poses a restriction on the possible strength of magnetic 
field in the X-ray emission region, 

.B < 0.1 [-Rpwn/10 13 cm] 1 G (16) 

In the model fit of Fig. 01 we have chosen the magnetic 
field strength B = 0.03 G and assumed the size of X-ray 
/ radio emission region 5pwn ^ 3 x 10 cm. 



5 Alternative mechanisms of TeV 7-ray 
emission. 

5.1 Bremsstrahlung. 

The IC model for the keV-to-TeV spectrum has a diffi- 
culty to explain the observed correlation of the radio, 
X-ray and TeV emission because the cooling and es- 
cape times of electrons emitting IC radiation in X-ray 
and TeV bands are different. Since the non-pulsed radio 
emission from the system is most probably related to the 
passage of the pulsar through the disk of Be star, an ex- 
planation of the observed correlation requires a physical 
mechanism which would explain the increase of the TeV 
flux during the disk passage. At least two mechanisms 
of interaction of the pulsar wind with the Be star disk 
can lead to the increase of TeV emission: bremsstrahlung 
and proton-proton interactions. 

As it is discussed above, the bremsstrahlung cool- 
ing time in the dense Be star disk © can be compa- 
rable to the IC cooling time both for the highest en- 
ergy electrons above TeV and for electrons with energies 
below 1 GeV (see Eqs. ©, ©)■ The bremsstrahlung 
cooling time can be comparable to the escape time from 
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Fig. 5 Comparison of bremsstrahlung and IC contributions 
to the 7-ray flux. See text for the values of parameters used 
for the model fits. 

the compact equatorial disk so that up to 10% of the 
power of the pulsar wind can be emitted in the form 
of bremsstrahlung radiation. Fig. [SJshows the fit for the 
7-ray spectrum of the system with a combination of IC 
and bremsstrahlung emission. The electron spectrum is 
supposed to be a cut-off powerlaw with the spectral in- 
dex p e — 2.5 and cut-off energy E cut — 20 TeV. Note 
that the EGRET upper limit imposes a restriction on 
the spectrum of electrons because the bremsstrahlung 
spectrum has the photon index r p h — p e ■ Assuming that 
the electron spectrum continues to lower energies with- 
out a break would violate the EGRET bound on the 
flux. A break at the energy E ~ 350 MeV was assumed 
in the electron spectrum in the model fit of Fig. [S] The 
break at this particular energy is naturally expected in 
the bremsstrahlung scenario, because below this energy 
the ionization loss dominates over the bremsstrahlung 
loss, which leads to the hardening of the electron spec- 
trum by Ap e ~ 1 at low energies. 



5.2 pp interactions. 

If the pulsar wind is proton-loaded, interactions of the 
pulsar wind protons with the protons from the dense 
Be star disk provide an additional source the TeV 7- 
ray emission. As is is discussed above, the "7-ray effi- 
ciency" of pp interactions is a factor of several lower than 
that of bremsstrahlung, but the relative contributions of 
bremsstrahlung and pp interactions into the 7-ray emis- 
sion depend on the proton-to-electron ratio of the PW. 
If the pulsar wind is mostly proton loaded, the ir° decay 
emission can dominate over the 7-ray emission from the 
pulsar wind electrons. An example of the fit to the TeV 
7-rayspectrum within the pp model is shown in Fig. 
We have assumed a powerlaw spectrum of protons with 
the spectral index p p = 2.6 for the model fit. Similarly 
to the case of bremsstrahlung, the EGRET upper limit 
on the flux imposes a restriction on the spectrum of the 
protons at several GeV energies. However, contrary to 
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Fig. 6 Broad band spectrum in the pp model. 7-ray emission 
is produced via neutral pion decays (red line) while radio and 
X-ray emission are, respectively, synchrotron and IC radia- 
tion from electrons/positrons produced in the charged pion 
decays. 



the bremsstrahlung case, the spectrum of ir° decays 7- 
rayemission has a low energy cut off below GeV energy 
and the restriction on the spectrum of the protons is 
rather weak. In fact, in the model fit of Fig.Elthe proton 
spectrum is described by a single powerlaw from GeV to 
TeV energies. It is important to note that in this case 
the luminosity of the source in the 10 GeV energy band 
is higher than in the TeV band. 

If bremsstrahlung or pion decay emission dominate 
the TeV flux, the TeV luminosity is just proportional 
to the density of the Be star disk at the location of the 
pulsar. This means that in this case the parameters of 
the Gaussian that approximately fits the TeV lightcurve 
(Fig. 01 should be identified with the equatorial plane 
position (6q) and width (A9q) of the disk. Within such 
interpretation the disk appears to be "thick" and the 
calculation of the phases of disappearance/reappearance 
of the pulsed emission has to be done by calculating the 
column density of the disk along the line of sight for dif- 
ferent positions of the pulsar. Fixing the phases of dis- 
appearance/reappearance of the pulsed emission to the 
known values, one can obtain a constraint on the radial 
density profile of the disk and/or inclination of the disk 
w.r.t. the line of sight. However, if inverse Compton emis- 
sion gives a contribution comparable to the pion decay 
or bremsstrahlung, the derivation of the disk parameters 
from the shape of the TeV lightcurve is not possible. 

The most important feature of the model with pp 
interactions is that it can explain the broad band spec- 
trum of the system from radio to TeV energy band. The 
idea is that the synchrotron and IC emission from the 
secondary electrons produced in the decays on charged 
pions is emitted in the radio and X-ray bands, respec- 
tively, while the bremsstrahlung emission from the sec- 
ondary electrons, emitted in the 7-ray band gives a sub- 
dominant contribution, compared to the tt° decay emis- 
sion. Fig.EJshows an example of the fit to the broad band 
spectrum of the system within the pp interactions model. 



Since the emission in radio, X-ray and TeV bands is pro- 
duced via one and the same process (pp interactions), 
the observed correlation of the radio, X-ray and TeV 
flux is naturally explained. Besides, the observed hard- 
ening of the X-ray spectrum during a several-day period 
following the moment of the entrance of the pulsar to 
the disk of Be star is explained by the low energy cut-off 
at ~ 100 McV in the spectrum of secondary electrons. 
Such a cut off arises (a) because of the kinematics of 
the pion decays and (b) because of the efficient Coulomb 
cooling of electrons with energies below 100 MeV during 
the escape from the Be star disk. 

The pp interaction scenario is attractive because of 
one more reason: in this model the overall energy bal- 
ance of the system is evident. Indeed, in the "purely elec- 
tronic" models it is not clear why the system is "radia- 
tively inefficient" : the spin-down luminosity of the pul- 
sar is ~ 10 36 erg/s, but the bolometric luminosity is just 
L < 10 erg/s, which accounts for no more than one per- 
cent of the spin-down luminosity. To the contrary, within 
"protonic" model one has to assume that proton-loaded 
PW carries a significant fraction of the spin-down power. 
As it is explained above, in the pp model the efficiency of 
conversion of the power contained in the protons into the 
7-ray emission is several percents, (see Eq. I|ll|l) which 
explains the 7-rayluminosity L 1 ~ 10 erg/s. 

The pp model can be readily tested with the future 
observations of the system in the 10 GeV energy band 
with GLAST. Indeed, from Fig.Elone can see that in the 
pp model the EGRET upper limit on the flux at 10 GeV 
should be close to the actual level of the 7-ray flux from 
the system. This means that the detection of the sys- 
tem during the periastron passage with a more sensitive 
instrument, like GLAST should not be a problem. 
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